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A method  is  summarized  for  the  precise  computation  of 
the  velocity  of  sound  in  a real  gas  and  of  its  variations 
with  experimental  conditions.  Specifically,  two  problems 
are  discussed:  (a)  the  inclusion  of  the  effects  of  depar- 

tures of  the  equation  of  state  from  the  perfect  gas  law, 
and  (b)  the  computation  of  specific  heats  from  spectro- 
scopic data.  The  method  is  applied  to  air  and  data  for 
the  computation  of  the  velocity  in  air,  under  pressures 
between  720  and  820  mm.  Hg,  and  temperatures  between  0°  C 
and  30°  C,  are  presented.  These  results  are  summarized  in 
the  form  of  correction  factors  for  the  reduction  to  standard 
conditions  of  velocities  measured  at  various  frequencies 
and  in  this  range  of  pressure  and  temperature. 
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Computation  of  the  Velocity  of  Sound  In  Gases 

by 

Preston  W,  Smith,  Jr. 

I 

INTRODUCTION 

The  recognized  factors  which  can  influence  the  phase  velocity 
for  propagation  of  sound  in  gases  under  various  experimental  con- 
ditions are  numerous  to  a point  of  annoyance.  However,  when  the 
gas  is  uniform  and  unbounded  and  the  temperature,  pressure,  and 
frequency  are  not  extreme,  the  situation  is  considerably  simplified. 

The  classical  expression  for  the  phase  velocity  of  an  acous- 
tic signal  of  small  amplitude  in  an  unbounded  gas  is  the  well- 
known  expression,^* 

c * wyoo)*  , u) 

where  Y is  the  ratio  of  specific  heats,  PQ  is  the  equilibrium 
pressure,  and  PQ  Is  the  equilibrium  density  of  the  gas.  The 
assumptions  implicit  in  this  equation  are  that  the  gas  is  "ideal" 
(i.e.,  the  equation  of  state  is  the  perfect  gas  law)  and  that 
the  sound  disturbance  is  an  isentropic  (reversible  adiabatic) 
process.  It  is  usually  assumed  that  Y is  independent  of  the 
experimental  conditions  $ it  follows  that  the  velocity  should 
vary  only  as  the  square-root  of  the  absolute  temperature. 

In  this  memorandum,  we  shall  consider  the  problem  of  the 
theoretical  calculation  of  the  velocity  of  sound  with  a pre- 
cision which  invalidates  some  of  the  preceding  assumptions. 

Pirst,  a method  will  be  outlined  for  the  incorporation  of  the 
effe  Is  of  the  "imperfection"  of  the  gas  (i.e.,  the  departure 
of  its  experimentally  determined  equation  cf  state  from  the 
perfect  gas  law).  Secondly,  the  means  for  the  computation  of 
the  specific  heat  from  spectroscopic  data  will  be  presented. 

No  consideration  will  be  given  to  the  effects  upon  the  velocity 

♦All  numbered  references  appear  at  the  end  of  the  report. 
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of  the  radiation  and  conduction  of  heat  or  the  viscosity,  3ince 
these  phenomena  are  of  importance  only  at  the  extremes  of  the  fre- 
quency scale*  A review  of  these  effects  can  be  found  in  a recent 

p 

paper  by  Markham,  et  al*  The  results  of  computations  for  air 
will  be  presented- 

From  the  results  of  theoretical  calculations  of  this  nature, 
measurements  of  the  velocity  of  sound  made  under  different  ex- 
perimental conditions  can  be  corrected  for  these  known  effects 
and  accurately  compared  at  a common  reference  condition* 

II 

THEORY 

The  distribution  of  pressure  in  a sound  field  in  an  unbounded 
medium  can  be  shown  to  be  the  solution  of  the  wave  equation,^ 

V2p  + («2/V2)p  =••  0,  (2) 

where,  in  lieu  of  a specified  equation  of  state  frr  the  medium, 
the  symbol  V2  stands  for  the  differential  quotient, 

V2  = (3p/ap).  (3) 

In  these  equations,  p is  the  instantaneous  variational  pressure, 
p is  ths  instantaneous  density  of  the  medium,  and  *>  is  the 
angular  frequency**  For  generality,  V must  be  considered  to  be 
a complex  number  Cthe  complex  velocity  ,q£  sound)  * However,  the 
situations  to  be  considered  in  this  memorandum  are  limiting  or 
asymptotic  cases  in  which  V is  real}  in  that  case,  it  is  readily 
recognized  as  the  phase  velocity  or,  more  commonly,  the  velocity 
of  sound. 


♦The  assumptions  implicit  in  the  mathematical  derivation  are; 

a.  the  pressure  is  a scalar  point  function  of  the  space 
coordinates, 

b.  there  are  no  sources  or  sinks  of  mass  in  the  region, 
c*  the  external,  or  "body,"  forces  on  the  medium  are 

constant, 

d*  the  pressure  and  density  disturbances  are  very  small,  and 
e=  the  time  derivative  of  V2  is  zero. 


I 


* 

& 
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1.  Velocity  of  Sound  In  Imperfect  Gases 

The  differential  quotient  for  the  square  of  the  complex 
velocity  [Eq  (3)3  can  be  evaluated  in  terms  of  fundamental 
thermodynamic  variables  from  the  equation  cf  state  for  the  gas. 
If  the  perfect  gas  law  is  assumed,  one  reverts  to  Eq=  (l)o 

The  method  re  will  use  for  evaluating  V was  apparently 

4. 

first  developed  by  van  Itterbeek  and  Keesom  for  determination 

of  the  real  gas  equation  of  state  from  sound  velocity  me a sure - 

c 6 

ments$  subsequently,  both  Richards'  and  Hardy  employed  similar 
expressions  in  the  consideration  of  the  inverse  problem  --  that 
with  which  we  are  concerned  — of  prediction  of  the  variations 
of  sound  velocity  from  equation-of-s tate  data. 

The  treatment  proceeds  from  the  fundamental  definition, 
[Eq.  (3)3,  and  by  means  of  standard  thermodynamic  relationships 
valid  for  pure  substances*  arrives  at  the  expression 

y2  = (fp5i  + (4) 

where  P is  the  total  pressure, 

T is  the  absolute  temperature, 

Cv  is  the  specific  heat  at  constant  volume,  and 
M is  the  molecular  weight  of  the  medium,. 


However,  further  manipulation  is  necessary  since  the  spe- 
cific heat  of  a real  gas  is  not  constant^  its  derivative, 
known  from  thermodynamics,  must  be  integrated  to  obtain  the 
value  of  Cy  at  the  temperature  and  pressure  of  interest,,  The 
equation  then  reduces  to 


V2  * &£  t 


1 + 


.x/.f 

C 

-I2  + a 


» 


(5) 


* The  term  “pure  substance”  is  used  here  in  its  thermodynamic 
sense.  Keenan  in  Thermodynamics  (Wiley,  194-1)  defines  it  ass 
"A  system  which  is  homogeneous  in  composition  and  homogeneous 
and  invariable  in  chemical  aggregation,, . " (p.  18).  The  re- 
striction is  necessary  in  order  to  eliminate  the  complexities 
of  phase  changes  and  chemical  reactions. 


i 

t 

\ 

i 
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where  R is  the  universal  gas  constant, 

M Is  the  molecular  weight  of  t'  gas, 

Cvo/R  is  the  ratio  of  the  specxxic  heat  at  constant 
volume  to  the  universal  gas  constant  in  the  limit  as 
the  pressure  approaches  zero,  and 

f,  g,  and  A are  dimensionless  functions  of  the  equi- 
librium values  of  P and  T.*  *♦ 


The  defining  equations  for  f,  g,  and  A (from  which  they 
can  be  evaluated  if  an  equation  of  state  relating  P,  P , and 
T is  known)  arex 


g 


A 


P2R2 


c -c 

v 


vo  _ T 
“ R 


(6) 


The  quantities  f and  g differ  only  slightly  from  unity,  and  A 
from  zero,  by  an  amount  determined  by  the  imperfection  of  the 
gas?  for  a perfect  gas  [Eq=  (5)]  reduces  to  [Eq.  (1)],  i.e., 
V2  - O'  p . 


The  significance  of  the  individual  factors  of  Eq.  (5)  can 
be  shown  to  be  as  follows? 


♦In  his  paper,  Hardy  (ref.  6.  above)  used  a function,  h,  re- 
lated to  our  A by  the  equations  (Cv/R)  - h(Cvo/R)  = (Cvo/R)  + A. 
The  difference  lies  primarily  in  the  fact  that,  whereas  A can 
be  evaluated  from  the  equation  of  state  alone, 'h  involves  also 
the  value  of  Cyo-  However,  Cvo  varies  with  temperature  inde- 
pendently of  A?  therefore,  the  present  notation  is  simpler  and 
preferable® 

♦♦The  factors  f,  g,  A,  and  Cyp/R  in  this  equation  are  all  dimen- 
sionless? any  consistent  set  of  units  can  be  used  in  the  compu- 
tations leading  to  their  evaluation.  The  units  of  the  resultant 
figure  for  the  velocity  of  sound  depend  wholly  on  the  factor, 
RT/Mc  If  it  is  desired  that  the  velocity  have  the  units, 
(cm/sec),  it  is  appropriate  to  express  T in  degrees  Kelvin, 

H in  grams  per  moles  and  R in  ergs  per  mole  « degree. 
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(a)  (RTf  TO  is  the  real  gas  isothermal  bulk  modulus  of 
elasticity?  -vCdP/dv)^,  where  v is  the  specific  volume, 

(b)  (g/f)  is  the  difference  of  the  specific  heats 
divided  by  the  gas  constant?  (Cp  - Cv)/R,  where  Cp  is  the 
specific  heat  at  constant  pressure. 

(c)  A is  the  increment  in  the  ratio,  Gy/R,  from  zero 
concentration  of  the  gas  to  the  experimental  concentration. 
Therefore,  the  quantity,  (Cvo/R)  + A,  is  the  actual  value  of 
Cv/R  as  defined  by  the  expressions  (3U/3T)V/R,  where  U is  the 
internal  energy  of  the  medium. 


(d)  The  whole  factor  in  square  brackets  in  Eq.  (5)  is 
therefore  the  real  y,  i.e,,  the  ratio  of  the  specific  heats 
evaluated  at  the  particular  temperature  and  concentration  for 
which  f,  g,  Cvo/R,  and  A are  computed. 

Of  course,  any  equation’ of  state  which  is  an  accurate 
analytic  representation  of  the  P,  p,  T interdependence  can 
be  used  to  compute  these  functions,  f,  g,  and  As  the  Beattie- 
Bridgman  equation^  is  one  of  the  best  available  for  this  pur- 
pose for  substances  in  vaporous  phase.  That  equation  expresses 
P as  a power  series  in  p,  the  coefficients  being  functions  of 
T$  it  involves  five  constants  determined  for  each  gas  by  a 
process  of  fitting  the  various  derivatives  of  the  pressure  to 
a best  match  with  data  measured  by  various  experimenters,  ^n 
its  general  form,  the  equation  can  be  written 


RdB. 


F = (RT)p  ♦ (RTB  - - -S2)p2  + (A  a - RTB  b - ) p3 


♦ (• 


o 

RdB  b 


'")  ■ 


v ( , 


*)  P 


where,  In  the  units  usually  employed, 

P is  the  absolute  pressure  in  atmospheres, 

R is  the  universal  gas  constant  (atmosphere  liter/mole 
degree) . 
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T Is  the  absolute  temperature  (degrees  Kelvin) , 
p is  the  molar  density  (moles/liter) , and 
A&  is  a constant  with  units  of  (atmosphere  liter^/mole^) , 
B , ni  are  constants  with  units  of  (liter/mole),  and 
d is  a constant  with  units  of  (degree-5  liter/mole ). 

This  equation  is  sometimes  conveniently  abbreviated  in  the 
/lrlal  form 


P “ ap  + (3p^  + YP8  4 Sp^  (8) 

where  the  four  coefficients  (a,  3,  Y>  5)  are  functions  of 
temperature  and  independent  of  density. 

Expressions  for  the  functions  needed  in  Eq.  (5)  for  the 

A, 

velocity  are  readily  evaluated  from  these  equations  as  power 
series  in  p.  In  the  general  case,, 


2dB 


2dB  b 


1 0 CUB.U 

g*  » 1 + (B„  4 *%)p  4 (-Bnb  4 — ^p2  +( £-)p3  , (9) 


S'L 


B b 0 

i 4 — 2 q _ _2_ 

-1-  r\  1*  1 ►* 


Jhe  following  approximate  formulas,  retaining  only  the 
first  non-"idealu  terras,  can  b<=  verified  from  the  expressions 
given  above, 

P = |f  Cl  4 (~B0  + b¥  + It3  (8a) 


f - 1 - 2(-Bq  + gj*  + ^)p 


2 An  ,A 

g/f  - 1 4 (-=8  ♦ 


RT 


A - p 
T3 


(9a) 
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In  sense  check  calculations  for  oxygen  at  room  temperature  and 
atmospheric  pressure,  the  results  of  the  approximate  formulas 
differed  but  slightly  from  the  results  of  the  more  accurate 
expressions.  The  errors  weres  less  than  *1  part  per  million 
in  p , +*  parts  per  thousand  in  (1-f),  -1  part  per  thousand  in 
(g-i),  and  -1  part  per  thousand  in  a. 

2c  Evaluation  of  the  Specific  Heat 

In  order  to  compute  the  velocity  of  sound  from  Eq„  (5) 
it  is  also  necessary  to  obtain  an  accurate  value  for  the 
specific  heat  at  infinite  dilution,  CyQc  It  has  frequently 
been  pointed  out  that  theoretical  computations  by  the  meth- 
ods of  statistical  thermodynamics,  with  data  obtained  by 
spectroscopic  analysis  of  molecular  spectra,  are  more  accurate 
than  direct  measurements.  However,  one  must  first  consider 
the  possible  differences  between  the  specific  heat  effective 
at  a given  sound  frequency  and  that  which  would  obtain  in  a 
quasi-static  process. 

Yazj&Slsa . yijthjre  opfinsj 

The  definition  of  the  specific  heat  can  be  written  as 

Cy  = (au/dT)v  (10) 

where  U is  the  internal  energy  of  the  medium,  T is  the  abso- 
lute temperature,  and  the  subscript  v signifies  a constant 
volume  process.  This  process  must  be  envisioned  as  a periodic 
variation  at  the  frequency  of  the  sound  wave.  Now  the  internal 
energy  consists  of  both  translational  energy  and  energy  associ- 
ated with  the  various  "excited”  (rotational  and  vibrational) 
states  of  the  molecules.  In  an  equilibrium  condition  a cer- 
tain small  fraction  of  the  molecules  exists  in  each  excited 
state,  the  proportions  being  dependent  on  temperature.  The 
periodic  variational  temperature  resulting  from  a sound  wave 
will,  therefore,  cause  a periodic  variation  of  the  proportions 
of  the  excited  molecules^  i.e.'  a periodic  variation  In  the 
internal  energy  associated  with  the  excited  states.  This 
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energy  is  in  addition  to  the  translational  energy  and  is 
therefore  seen  to  modify  the  effective  specific  heat.  How- 
ever, the  adjustment  in  the  proportions  of  molecules  in  the 
various  states  is  not  instantaneous^  it  can  be  shown  that 
there  is  a "time  constant"  associated  with  the  response 
mechanism  which  introduces  a phase  difference  between  the 
temperature  variations  and  the  resultant  energy  variation. 
Therefore  the  effective  specific  heat,  Cy,  and  the  velocity, 

V,  will  each  be  complex.  In  two  cases,  Cy  and  V will  be 
reals  first,  wnen  the  sound  frequency  is  very  low  and  th9 
gas  can  be  considered  to  be  moving  through  equilibrium  states 
only$  secondly,  when  the  sound  frequency  is  so  very  high 
that  there  is  no  time  for  an  adjustment  in  the  proportions 
of  molecules  in  the  various  states.  Between  these  two  limit- 
ing cases,  the  effective  specific  heat  will  have  decreased  by 
an  amount  equal  to  the  contribution  of  the  energy  in  the 
excited  states  to  the  quasi-static  specific  heats  the  velocity 
will  have  increased  correspondingly. 

In  the  general  case  it  is  found  that  the  transition  for 
all  the  states  in  any  one  "degree  of  freedom"  occurs  at  a 
single  characteristic  frequency,  although  these  frequencies 
may  be  (and  usually  are)  different  for  different  degrees  of 

O 

freedom.  The  picture  that  develops  is  the  following.  The 
quasi-static  specific  heat  (that  effective  at  zero  frequency) 
includes  terms  attributable  to  each  of  the  translational, 
vibrational,  rotational,  and  electronic  degrees  of  freedom 
of  the  molecules.  As  the  frequency  increases,  these  terms 
successively  drop  out  from  the  effective  specific  heat.  Thus 
the  velocity  should  increase  in  a sequence  of  steps.  A careful 
analysis  shows  that  the  velocity  makes  a smooth  transition 
between  these  steps,  but  at  frequencies  far  from  the  transition 
region  t}ie  value  of  the  velocity  is  essentially  that  of  the 
"step."  (See  Fig.  (1).)  We  shall  be  concerned  here  with  the 
computation  of  the  values  of  these  steps.  Therefore  we  wish  to 
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evaluate  the  contribution  of  each  of  the  degrees  of  lreedom 
to  the  specific  heat*  More  specifically,  only  the  limiting 
values,  as  the  pressure  approaches  zero,  must  be  known  (see 
Eq*  (5))^  therefore,  the  gas  can  be  considered  an  ideal  gas 


I 

fre quency 


Fig.  L Variation  of  Velocity  with  Frequency* 

in  this  section*  The  determination  of  the  transition  frequen- 
cies for  the  varJous  degrees  of  freedom  is  an  entirely  distinct 
problem  which  is  not  considered  in  this  memorandum*  It  is,  at 
present,  approached  most  readily  by  the  experimental  measure- 
ment of  the  dependence  of  the  velocity  and  attenuation  of 

O 

sound  upon  frequency* 

Computation  from  Spectroscopic  Data 

There  are  a number  of  papers  and  books  which  discuss,  in 
varying  degrees  of  completeness  and  complexity,  the  computation 
of  specific  heats  of  gases  from  spectroscopic  data*  however,  it 
is  appropriate  to  summarize  here  the  basic  theory  and  the  meth- 
ods used  in  this  memorandum*  Since  this  report  is  primarily 
concerned  with  the  gaseous  mixture,  air,  attention  will  be 
centered  on  its  component  gases  which  are  mostly  dj atomic*  The 
many  complications  of  detail  which  obtain  for  polyatomic  gases 

Q 

and  at  extremes  of  temperature  are  thereby  avoided*  7 

The  fundamental  theoretical  expression  for  the  specific 
heat  at  constant  volume  of  a pure  gas  is1-0 
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Cy/R  = T'2  32  -ea  f(T)/3(l/T)2,  * (11) 

where  f(T)  is  the  partition  function.  The  partition  function 
for  a gas  depends  on  the  temperature  and  the  energy  levels  of 
the  various  states  of  the  molecules,,  Its  analytic  expression 
is9  ^ — 7 

f(T)  = AT3/2  Pi  e ^ ^ (12) 

i 

where  A is  a constant  independent  of  temperature,  p.  is  the 
relative  statistical  weight  of  the  itn  state,  and  E4  is  the 
energy  (in  units  of  degrees)**  of  the  i state,  The  statis- 
tical weight  is  the  number  of  distinct  configurations  of  the 
molecule  which  have  the  same  energy,  E^ , and  therefore  need 
not  be  distinguished  in  the  sum  for  the  partition  function. 

It  is  readily  verified  from  Eqs.  (11)  and  (12)  that  neither 
the  constant  multiplier  A nor  a change  in  the  datum  for 
energy  (the  same  change  for  all  states)  will  affect  the  value 
of  the  specific  heat. 

If  the  equation  (12)  is  substituted  into  Eq,  (11),  the 
results  can  be  written: 


where 


Cy^H 

.trans/g 


C*nt/R 


.int 


« (C*rans/R)  + (C*nt/R),  (13) 

= T"2  d2  -6i(T3/2)/d(l/T)2  = 3/2, 

= t"2  a2  fln'c/a(i/T)2, 


♦The  subscript  "o".  indicating  evaluation  at  zero  concentration, 
has  been  omitted  rrom  Cy  throughout  this  section, 

♦♦The  energy,  in  degrees,  of  a particular  state  can  be  obtained 
from  the  corresponding  wave  number  (in  car^)  by  multiplication 
by  the  constant,  hc/k  where  h is  Planck's  constant,  c is  the 
speed  of  light,  and  k is  Boltzmann's  constant.  The  two  units 
are  entirely  equivalent. 


i 
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The  first  term,  ecual  to  3/2,  is  the  contribution  of  the  trans- 
lational energy  z the  specific  heat.  The  second  term  measures 
the  contribution  of  the  ’•internal"  degrees  of  freedom  or,  in  the 
terminology  used  earlier  * **  the  excited  states  of  the  molecule. 

The  evaluation  of  this  term  is  the  major  concern  of  this  section. 

It  is  seen  that  the  internal  partition  function t f 1 * , is  formally 
determined  by  the  energy  levels,  E^,  and  their  corresponding  weights, 
p^«  These  energy  levels  (in  units  of  the  wave  number,  cm-1)  can 
be  directly  measured  spectroscop Lcally  from  the  absorption  and 
amission  spectra  of  the  gas.  It  is  essential  to  the  systematic 
computation  of  the  specific  heat,  however,  that  one  employ  the 
theoretical  forms  for  the  differences  between  successive  energy 
levels  as  determined  by  quantum  mechanics  with  values  of  the 
necessary  constants  determined  by  a best  match  with  the  experi- 
mental spectroscopic  data. 

In  the  usual  convention,  the  total  energy  of  a molecule  in  a 
given  s i S expanded  in  the  form 

E = Ee(e)  + Ev(e,v)  + B^(e,v,J)  , (14) 

where  E (e)  is  the  energy  level  of  the  e-^  electronic  state  of  a 
molecule,  in  the  hypothetical  case  where  there  is  no  vibration  or 
rotation,  E_(e,v)  is  the  additional  energy  which  the  molecule 

7 w 7 .yy.  ^y. 

would  have  if  it  were  in  chc  v — vibrational  level  of  the  e^ 
electronic  state,  without  rotation,  and  E (e,v,j)  is  the  energy 
possessed  by  a molecule  in  the  rotational  level  in  excess  of 
the  total  energy  without  rotation.  Furthermore,  the  statistical 
weight,  p1?  can  be  factored  into  the  pruduct  of  a weight  for  the 
rotational  level,  Pr(j),  and  a weight  for  the  electronic  state, 

P«(e) . 

w 

The  internal  partition  function  can  now  be  written 

iL  e 

e=0  v-0 


,int  _ 


■2- 


_E  (e) 
((e)e  — ijr~ 


-V'e,v) 

T 


.Er(«,T,. 


(15) 
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In  many  cases,  in  gases,  the  couplings  between  the  rotational, 
vibrational,  and  electronic  states  are  sufficiently  small  that 
their  effects  on  the  partition  function  can  be  neglected . Mathe- 
matically, this  statement  means  that  if  the  values  of  the  sums 

X and  Z are  relatively  independent  of  ihe  other  quantum 
j=0  v=b 

numbers,*  their  dependence  on  those  numbers  can  be  neglected. 

This  is  particularly  true  when  the  value  of  Ev  (or  E^)  increases 
greatly  with  the^uantum  number  so  that  the  multiplicative  factor, 

e (or  e 8 ) decreases  sharply. 

When  couplings  between  degrees  of  freedom  are  negligible, 
their  individual  contributions  to  the  specific  heat  can  be  shown 
[from  Eqs.  (13)  and  (15)3  to  be  independent  and  additive.  Thus 
we  can  write 

Cvint/R  = (Cvel/R)  + (Cvvib/R)  + (Cvrot/R)  , (16) 

when  the  individual  terms  are  to  be  evaluated  bv  means  of  Eq.  (11) 
from  the  partition  functions  for  the  individual  degrees  of  freedoms 

„ _ Ee(e)-Ee(0) 

fel  c P*(°)  + Pft(«)  E * 

9 e*l  6 

Ev(v)-Ev(0) 

e T (17) 

ErO) 

tr^  = PrU)  e “ 1 . 

j=0  r 


f 


,vlb  _ 


= 1 + 


z 

v=l 


V 


♦The  identifying  index  for  the  electronic  state,  e,  is  not  prop- 
erly called  a quantum  number^  that  term  is  being  misused  here 
in  order  to  simplify  the  discussion.  A further  complication 
of  nomenclature  arises  in  the  case  of  polyatomic  molecules 5 
the  sum  over  v must  then  be  replaced  by  a double  sum  since  two 
numbers  are  necessary  to  identify  a single  vibrational  state? 
a modal  index  and  the  quantum  number  identifying  the  level 
within  the  mode. 
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Method s of  computation , 

There  are  two  particularly  useful  methods  for  evaluating  the 
contributions  to  the  specific  heat  from  these  series.  The  value 
of  (C  vib/R)  can  be  reduced  tc  the  closed  form11 


Cyvlb/R  = x2/sinh2x  9 

where  (18) 

x * [ E^.  ( 1 ) -Ev ( 0 ) 3 /2T  . 


This  expression  is  derived  under  the  fairly  accurate  assumption 
trial  the  suc^os.iiva  vibrational  energy  levels  are  evenly  spaced 
(the  "harmonic  vibrator"  assumption).  However,  the  expression 
is  also  a valid  approximation,  without  the  harmonic  assumption, 
i£  the  temperature  is  low  enough  that  only  the  lowest  two  vib- 
rational levels  are  significantly  populated  (i.e.,  if  e ~*x  is 
negligible).  Therefore,  Eq.  (18)  is  very  accurate  at  room  tem- 
pervitur-e  but  is  somewhat  in  error  at  high  temperatures.* 

The  most  accurate  method  for  evaluating  any  of  the  contri- 
butions to  (Cv/R)  is  computation  directly  from  the  series  (Eqs. 
(17))  by  means  of  Eq.  (11).  The  labor  is  readily  adapted  to 
machine  calculation.  It  can  be  verified  that,  if  the  partition 
function  has  the  form 


the  contribution  tc  the  specific  heat 


VR  = 


f f"  -•  (f’)2 

72~ 


where 


f'  P1(Ei/T)e  1 


(19) 


♦This  equation  is  also  valid  for  the  contributions  of  each 
of  the  vibrational  modes  of  a p:xyatomic  molecule,  if  the 
modes  can  be  considered  "harmonic  vibrators." 
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f"  Pi  (Bj/Tr  s 

This  method  of  calculation  must  be  used  to  evaluate  the  contri- 
bution of  the  rotational  degrees  of  freedom  and  can  be  used  for 
each  of  the  internal  degrees  of  freedom. 

Formulas  for  the  enerev  levels  and  statistical  weights 

There  remains  only  the  problem  of  determining  the  values  of 
the  energy  levels  and  the  statistical  weights  for  the  various 
degrees  of  freedom.  As  discussed  earlier,  ws  shall  use  the  theo- 
retical formulas  determined  by  quantum  mechanics. 

The  electronic  contribution  to  the  specific  heat  is  vanishing 
small,  in  the  gases  of  interest,  at  all  reasonable  temperatures 
because  of  the  large  values  of  the  differences  in  energy* levels, 

(E  (e)~E  (0) ) o [See  Eq,  (17).]  In  oxygen,  in  which  the  effect 
is  larger  than  in  most  gases,  the  first  electronic  state,  XA, 
above  the  ground  state,  , contributes  less  than  1 part  in 
1C4  for  temperatures  less  than  SGO0?:.1^ 

The  spacing  of  the  energy  levels  of  the  vibrational  degree 
of  freedom  can  be  written  approximately  as 

Ev(v)  - Ev(0)  = (we  - V>e)v  ' (xea)e)  ^ » (20) 

where  o>e  is  a wave  number  (reduced  to  degrees),  and  xq  9s  a small, 
positive  constant  included  to  compensate  for  the  anharmonicity  of 
the  seccessive  energy  levels.^4  The  values  of  and  xQ  are 
tabulated  in  references  on  molecular  spectra. 


The  rotational  energy  levels,  for  a nonrigid  molecule  in  the 
v—  vibrational  level,  are  given  by  the  equation"1' 


Er(v,J)  = Bv  JU+D  + Dv  J2U+1>2 


(21) 


- Be  - a(v  ♦ §) 

Dv  = De  + e(T  * V 


where 


9 


Tlfi?9  - 1 5- 

Here,  tne  small  factor  Dy  corrects  lor  the  effect  of  the  non- 
rigidity  of  the  molecule  in  rotation,  this  effect  is  sometimes 
called  ’‘centrifugal  stretching  . " The  small  factors  a and  3 
correct  for  the  coupling  between  the  vibrational  motion  and  the 
rotation.  Were  this  formula  to  be  used  in  calculation,  the 
partition  function  for  the  vibration- rotation  levels  could  not 
be  separated  into  products  as  in  Eq.  (15)  nor  could  the  contri- 
butions of  vibration  and  rotation  to  the  specific  heat  be  separated 
as  in  Eq.  (16) « However,  since  the  coupling  is  small,  its  effect 
on  the  specific  heat  is  quite  negligible  at  all  but  extremely 
high  temperatures.  It  is  therefore  sufficient  to  use  Eq.  (21) 
with  v = 0.  The  quantities  usually  tabulated  arc  tne  constants 

B-  = -(a/2)  and  a . The  various  constants  are  related  to 

° e 

one  another  by  the  equations  ' 

= - 4 BeV< 

3/'De  1 (ce/24  Be)(a/Be)2  + 5 (a/B#)  • 8 I(  . (22) 

By  these  equations,  all  constants  can  be  determined  from  the 
tabulated  values  and  the  energy  levels  computed. 

It  is  still  necessary  to  determine  the  statistical  weights, 
d ( j ) , for  the  rotational  levels.  For  the  gases  of  major  interest 

■T  1 iZ 

in  this  report,  these  weights  can  be  written 

Pr(J)  * o j ( 2 j ♦ 1)  , (23) 

where  the  value  of  g^  is  different  for  odd  values  of  j frem  that 

for  even  j’s.  The  ratio  (g.  „)  is  a constant  for  the 

u onn  even 

gas.  All  the  equations  necessary  for  the  computation  of  the  rota- 
tional contribution  to  the  specific  heat  have  now  be^n  summarized.*  ** 

♦Strictly  speaking  Eqs.  (21)  and  (23)  apply  only  to  singlet  elec- 
tronic states  such  as  the^ground  state  of  In  the  ground 

state  of  O2,  the  triplet  the  energy  levels  and  statistical 
weights  vary  slightly  between  the  three  states.  Fortunately,  this 
distinction  is  not  significant  except  at  very  low  tempera tures^  at 
room  temperatures  and  higher,  the  ground  states  of  O2  can  be  con- 
sidered to  be  a single  degenerate  state  of  weight , 3 « 

**It  is  frequently  assumed  that,  at  room  temperature,  the  rotational 


TK29 


-16- 


Specific  Heat  of  Mixtures  of  Pure  Gases 

The  preceding  comments  have  been  concerned  with  the  cal- 
culation from  spectroscopic  data  of  the  specific  heat  of  a pure 
gas  at  zero  concentration,  Cv0./R.  Since  a gas  at  zero  concen- 
tration can  be  considered  to  be  an  ideal  gas,  the  specific  heat 
of  a mixture  can  be  computed  by  the  usual  formula, 

(Cv0/R)tot  ~ (Cvo/R)i  ’ 

♦“Vi 

where  x,  is  the  mole  fraction  (part  by  volume)  of  the  i 

1 17 

component- = 


states  of  the  molecule  are  so  fully  excited  that  the  rotational 
degrees  of  freedom  are  behaving  entirely  classically  and  that 
the  contribution  of  each  to  the  specific  heat,  Cv/R,  is  there- 
fore equal  to  1/2,  This  is  strictly  true  only  for  the  theoretical 
case  of  a rigid  rotator.  The  calculations  performed  for  this 
report  indicate  that  the  effect  of  centrifugal  stretching  in- 
creases the  value  of  CyroVR  for  oxygen  by  two  parts  per  thousand 
at  room  temperature.  The  effect  increases  with  increasing  tempera- 
ture. 
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III 

COMPUTATIONS  AND  RESULTS  FOR  AIR 
1.  Computations  for  Imperfections  of  Air 

The  equation  for  the  square  of  the  velocity  of  sound  in 
an  Imperfect  gas  was  der  if**  in  a previous  section.  [See  Eq.  (5), 
p,  4.3  The  effect  of  the  imperfection  of  the  gas  Is  contained 
in  three  functions:  f,  g,  ana  A which  must  be  evaluated  from  the 

equation  of  stuto.  Expressions  for  the  three  runetlons  have  been 
derived  [Eq.  (9),  p.  6]  in  terms  of  the  temperature,  the  density, 
and  the  five  constants  AQ,  a,  B , b,  and  d which  appear  In  the 
Beattie-Bridgman  equation  of  state. 

The  process  of  computation  of  these  functions  requires  same 
comments.  First  it  is  to  be  noted  that  the  independent  variables 
in  Eq.  (9),  as  in  the  equation  of  state  [Eq.  (7)3,  are  the  moLar 
density  and  the  absolute  temperature.  Since  the  measurable  quan- 
tities are  pressure  and  temperature,  it  is  necessary  first  to 
compute  the  molar  density,  p,  [by  a process  of  successive  approxi- 
mation, using  Eq.  (7)3  for  se‘ected  values  of  P and  T.  The  unit 
of  pressure  is  the  normal  atmosphere  so  that  conversion  is  necessary 
when  the  measurement  is  made  in  millimeters  of  mercury.  Finally 
the  values  of  the  universal  gas  constant,  R,  and  the  ice  point, 
T0(°K)used  by  Beattie  and  Bridgman  in  the  computation  of  the 
constants  Aq,  a,  Bc,  b,  and  d are  R = 0cC8206  atmosphere  liter/ 
mole  degree  and  TQ  = 273-13  degrees  Kelvin,  respec tively;  these 
figures  must,  of  course,  be  employed  in  all  further  calculation? 
with  these  five  constants,  namely,  in  the  determination  of  the 
values  of  f,  g,  and  A.  However,  in  Eq.  (5)  for  the  square  of 
the  velocity  of  sound,  the  absolute  temperature  appears  explicitly; 
here  the  best  contemporary  value  for  the  ice  point  should  be  used. 

1 P 

The  value  TQ  = 273.16  degrees  Kelvin  has  been  used  in  this  report. 
The  values  of  the  desired  functions,  f,  g,  and  A,  and  the  ratio 
g./f,  can  now  be  calculated  from  Eq.  (9)  for  any  particular  gas  or 
vapor  for  which  the  five  constants,  Aq9  a,  B , b,  d,  have  been 
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deternunedr  Mixtures  cf  these  gases  can  also  be  treated  oy  a 

iq 

method  developed  by  Beattie, 

Table  I presents  the  result  of  such  calculations  of  the 
functions  f,  g/'f,  and  A,  for  air  in  the  region:  0 to  30°C, 

720  to  820  mm.  Hg-  the  values  used  for  the  constants  of  the 
Beattie-Bridgnan  equation  were 

Aq  * 1,3012  Bq  = 0.04611  d = 4.34  X 104 

a = 0,01931  b = -0.01101 

The  results  contained  in  Table  I,  with  the  knowledge  of 
the  specific  heat,  Cv0/R,  are  sufficient  to  compute  the  ratio 
of  velocities  of  sound  at  different  conditions  of  temperature 
and  pressure.  However,  in  order  to  compute  an  absolute  value 
for  sound  velocity  from  these  theoretical  considerations,  it 
is  necessary  as  well  [see  Eq.  (5)3  to  know  the  value  of  the 
molecular  weight,  IS.  M is  not  determinable  from  the  equation 
of  state  alone 5*  indeed,  it  never  appears  therein  when  the  volume 
unit  is  molar  (e„  g, , liter/mole) . It  can,  however,  be  determined 
simply  from  molar  density,  p (mole/liter),  and  the  density  proper, 
p'  (graras/liter) , if  both  are  known  at  the  same  experimental 
conditions,  for,  by  definition,  p'  =p  M.  Furthermore,  the  value 
of  the  molar  density  is  readily  computed  from  the  equation  of 
state  [Eq.  (7)3=  Therefore  the  problem  reduces  to  determining 
the  density  in  grams/liter. 

*1  can  find  no  satisfactory  answer  in  the  literature  to  the 
problem  of  how,  from  experimental  data  in  which  the  unit  of 
volume  was  Dresumably  in  cgs  units  (cnw,  liter,  or  specific 
volume  units  of  cm3/gm,  liter/gm,  or  equivalents),  the  con- 
stants for  Eq.  (7)  could  have  been  determined  in  molar  units 
without  simultaneously  determining  the  value  either  of  the 
molecular  weight,  SJ,  or  of  the  individual  gas  constant  which  is 
the  ratio  R/M.  Indeed,  in  their  first  paper,  £.  Am.  Chem.  Soc . 
49.  1665  (1927),  Beattie  and  Bridgman  present  resuTts  for 
ethyl  ether  and  carbon  dioxide  which  give  both  the  values  of 
the  five  constants  and  the  value  of  the  ratio,  R/K;  none  cf 
the  other  papers  do  so. 
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Unfortunately,  nc  single  value  for  the  density  of  air  can 
be  stated  for  the  simple  reason  that  the  word  "air"  comprehends 
a whole  class  of  vaporous  mixtures  cf  sensibly  different  composi- 
tions. One  can  even  eliminate  the  well-known  variations  in  the 
proportions  of  water  vapor  and  carbon  dioxide  and  in  the  relative 
proportions  of  Op  and  Np  at  different  altitudes;  there  remains 
a fluctuation  of  density,  of  approximately  *0.0004  gram/liter, 

more  or  less  dependent  on  the  barometric  pressure  at  the  time 

21 

the  sample  was  obtained  for  measurement.  This  variation  is 
ascribed  to  changes  in  the  relative  proportions  of  oxygen  and 
nitrogen  brought  about  by  mass  movements  of  gas  within  the  at- 
mosphere — an  atmospheric  "stirring”  in  times  of  changing 
weather  conditions.  (This  amount  of  variation  in  density  causes 
a corresponding  variation  in  the  velocity  of  sound  of  about 
*0.05  m/s.)  Perhaps  as  satisfactory  an  answer  as  possible  is 
that  the  density  of  dry  air  with  0.03  per  cent  C0p,  at  sea  level, 
should  average  (over  months)  about  1.2930  X 10“^  grams/cm*  under 
standard  conditions  of  temperature  and  pressure.  ’ J Assuming 
this  value  to  be  exact,  the  molecular  weight  of  air,  computed 
by  the  method  described  above,  is  M = 28.9*2 7 grams. 

It  must  be  mentioned,  however,  that  these  slight  variations 
in  composition  of  "air"  and  in  the  molecular  weight  of  the  mix- 
ture should  not  significantly  affect  the  values  of  the  functions 
f,  6/'f,  and  A used  to  correct  measurements  to  standard  conditions 
of  temperature  and  pressure. 

2.  Computations  for  Specific  Heat  of  Air 

A summary  was  given  in  the  last  chapter  of  the  methods 
and  formulas  used  to  compute  the  specific  heats  of  pure  gases 
and  mixtures.  We  shall  present  here  the  values  of  the  experimental 
constants  and  the  results  of  the  computation. 

For  the  purposes  of  this  section,  ciT'  is  considered  to  be 
a mixture  of  four  gases  mixed  in  the  following. ..proportions  by 
volume : 
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N2,  78o 09^5  02,  20.95#*  A,  0.93#*  C02,  0.03#. 

Although  sore  variation  in  the  values  of  the  proportions  of 

these  gases  aonears  in  the  literature,  these  figures  are  among 

24 

the  most  recent . 

Two  constants  are  necessary  to  evaluate  the  exponents, 

E/T,  from  the  spectroscopic  data  given  in  (cm)"^  units  * the 
values  used  are: 

absolute  temperature  of  the  ice  roint;- 

Tq  = ?73.16°K; 

energy  level  cnnvprston  factor  (cm  x to  deg): 

= I.438337  deg°cm. 

The  data  used  for  the  individual  gases  follow: 

Argon  is  taken  as  a monatomic  gas  with  energy  of  translation 
only 5 Cv/R  = 1.9000. 

Carbon  dioxide  is  taken  as  a linear  polyatomic  molecule  with 

full*  excited  rotational  states  and  only  slightly  excited  higher 

vibrational  levels.  Precisely,  account  is  taken  of  the  energy 

of  the  first  vibrational  levels  above  the  ground  state  in  each 

of  four  modes  for  which  the  energy  levels  are  taken  as: 

Ev(*,v)-Ev(0,q)  = 9*54,  9.54,  18.9,  33=6  degrees, 

respectively.  Equation  (17)  ia  then  used  to  evaluate  the  cen- 
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tribution  of  each  mode  to  the  specific  heat.  ' The  result  need 
not  be  too  accurate  since  the  mole  fraction  of  C0o  in  air  is 
so  small. 

Oxygen  is  treated  as  a diatomic  molecule  in  a degenerate 
electronic  state  (^51).  The  contribution  of  the  next  highest 
electronic  state  ) to  the  specific  heat  is  neglected.  The 
small  coupling  between  vibrational  and  rotational  states  is 
ignored  and  the  vibrational  component  of  the  specific  heat  was 
computed  at  5“degree  intervals  from  0°C  to  3C°C  by  the  use  of 
Eq  - (18).  These  data  were  checked  by  direct  computation  from 
t he  series  expansion  for  the  partition  function  at  temperatures 
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of  0°C  and  20°C$  no  significant  difference  was  detected.  The 

rotational  specific  heat,  C r0^/R,  was  computed  at  0°C9  10°C, 

and  20UC  from  the  appropriate  series  by  Giauque;s  method  [Eq.  (19) 3 • 

The  values  of  the  spectroscopic  constants  used  for  these  com- 
26 

putations  are: 
vibrational  contribution: 


(in 


electronic  state) 


coe  = 1564.91  cm  , 
rotational  contribution: 

Bo  ' ' ! ’ 1-U8  cm 

a = 0.016  cm 


a)  x = II  064n  cm” 
e e 

(in  electronic  state; 
pr  = ( 2 j+1 ) , J odd 
=0  , j even 


Computations  for  Nitrogen  were  made  in  a manner  similar 

to  that  for  oxygen.  The  vibrational  specific  heat,  Cvvi^/R, 

was  computed  at  intervals  of  5°C  from  0cC  to  30°C,  using  Eq.  (18)0 

The  rotational  specific  heat,  Cvr0^/R,  was  computed  at  0°C  and 

20°C  from  the  series  expansion  for  the  partition  function  by 

Giauque's  method.  The  values  of  the  parameters  used  in  these 

27 

compn tat ions  are;  ' 
vibrational  contribution: 

toe  = 2359 06O  cm  ^ 

rotational  contribution: 

Bo  = se  - f*  1.992  cm 

a = 0.018  cm' 

In  Table  II,  the  computed  valxies  cf  specific  heats  are 
presented  for  each  gas  and  for  the  mixture,  air.  The  specific 
heats  for  air  are  also  given,  omitting  the  vibrational  con- 
tributions of  various  combinations  of  the  component  gases 5 

these  would  be  the  effective  specific  heats  at  various  high 
oA  “ 

frequencies/-  The  few  values  for  snecific  heats  riven  by 
29 

Hardy  show  no  significant  difference  from  these  results. 


1*T7 

(in  electronic  state) 

co,x  = 14.445  cm-1 
e e 

(in  electronic  state) 
pr  = 2(2i+l),  1 odd 

= (2J+1)  J even 


TU29 


-22- 


3 o Computed  Velocities  for  Air 

The  results  of  the  computations  of  the  previous  two  sections 
can  now  be  combined  and  theoretical  values  computed  [by  Eq.  (5), 
p.  4)]  for  the  velocity  of  sound  m air  and  its  variations  with 
temperature , pressure  and  frequency. 

The  velocity  of  sound  In  air.  at  the  reference  conditions 
generally  accepted  as  standard,  is  thereby  computed  to  be 
331*45  ra/s.  These  reference  conditions  are:  temperature,  0°C; 

pressure,  760  mm  Hg;  proportion  of  CO^  by  volume,  0.03  & humidity 
zero;  frequency,  sufficiently  low  that  the  specific  heat  should 
assume  its  quasi-static  value. 

The  computed  variations  of  the  velocity  of  sound  with  temper 
ature,  pressure,  and  frequency  are  nresented  In  the  accompanying 
figures  (Fig.  2)  in  the  form  of  the  fractional  correction  which 
should  be  applied  to  experimental  measurements  to  reduce  them  to 
reference  conditions.  The  correction  factor  given  is  the  amount 
by  which  the  ratio  of  the  velocity  at  reference  conditions  to 
the  velocity  at  stated  conditions  exceeds  unity.  The  curves  are 
lines  of  constant  correction  factor  (and  therefore  of  constant 
velocity) 5 they  are  plotted  as  functions  of  temperature  and 
pressure.  Four  charts  are  miven  representing  different  fre- 
quency regions.  The  four  charts  in  Fig.  (2)  represent  data  com- 
puted with  effective  specific  heats  for  air  which  omit  various 
contributions  from  the  internal  degrees  of  freedom  of  the  com- 
ponent gases.  Thus  in  Fig.  (2a)  (’’Low  Frequencies"),  no  con- 
tribution is  omxtted;  in  Fig.  (2b)  ("Intermediate  Frequancies  I") 
the  contribution  of  the  vibrational  degrees  of  freedom  of  is 
omitted;  in  Fig.  (2c)  ("Intermediate  Frequencies  II"),  the  con- 
tributions of  the  vibrational  degrees  of  freedom  of  0~  and  N ^ 
are  omitted;  in  Fig.  (2d)  ("High  Frequencies")  all  vibration 
contributions  are  omitted  and  the  effective  specific  heat  in- 
cludes only  translational  and  rotational  terms. 
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The  frequency  regions  in  which  the  charts  are  valid  are 
those  far  from  the  transition  frequencies  where  the  velocity  is 
changing  most  rapidly,,*  [See  Fig.  (1),  p.  Unfortunately 

the  transition  frequencies  are  not  definitely  determined  for  all 
gases.  Knotzel w has  indicated  that  the  transition  frequency 
attributable  to  the  vibration.al  degree  of  freedom  of  02  is  rot 
above  about  40-60  c/s  in  dry  pure  oxygen  dry  air,  however  he 
was  unable  to  fix  its  position  definitely.  Measurements  of  the 
ratio  of  specific  heats  of  02  made  by  Koehler7^  at  a frequency 
of  about  1 c/s  tend  to  indicate  tnat  that  frequency  is  still  some- 
what above  the  transit ion,  unfortunately,  the  experimental  probable 
error  is  too  great  to  allow  much  dependence  upon  this  deduction. 
Other  measurements  have  not  given  significant  results  for  dry  air.** 
There  is  no  information  atout  the  transition  frequency  attributable 
to  the  vibrational  degree  of  freedom  of  Some  of  the  transition 

frequencies  for  the  three  vibrational  modes  of  CO2  have  been 
located.  The  deformation  mode  is  most  important  of  these  and  its 
transition  frequency  (in  dry,  pure  C02)  is  about  17  kc/s.^2 

The  charts  in  Fig.  (2)  have  been  computed  on  the  assumption 
that  the  transition  frequency  for  the  vibrational  degree  of 
freedom  of  lies  at  some  unknown  point,  "X",  above  that  for  02« 
Under  this  assumption,  Fig„  (2a)  is  applicable  below  the  transition 
frequency  for  C~i  which  is  not  higher  than  40-60  c/s;  Fig.  (2b) 


*The  numerical  value  of  the  term,  "far  from,"  depends  on  the 
desired  accuracy.  It  can  be  shown  from  the  standard  equa- 
tions for  the  phase  velocity  in  the  transition  region  [V.0. 

Knud  sen.  Acoust.  Soc.  Am.  £,  199  (19^5)1,  that  88  per  cent 
of  the  total  change  in  velocity  takes  place  in  a frequency 
band  of  two  octaves  on  either  side  of  the  trai  'ition  frequency. 
Such  a factor  of  four  in  frequency  is  therefore  usually  quite 
adequate . 

**0f  course,  this  transition  frequency  has  been  determined  in 
moist  air  and  moist  oxygen,  in  which  cases  it  is  much  higher. 
[See  V.  C.  Knud  sen-  jL>  Acoust.  goc . Aa*  199  ( 19?  5)3-  How- 
ever, neither  the  aata  on  the  effect  of  Imperfection,  of  the 
gas  nor  the  specific  heats  apply  directly  to  moist  air. 
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is  applicable  between  this  point  and  "X";;  Fig.  (2c)  is 
applicable  between  "X”  a~d  the  trarsition  frequency  of  C0o 
at  17  kc/sj  finally  Fig.  (2d)  is  applicable  in  the  frequency 
region  above  17  kc/s0  The  upper  lit.,  11  of  applicability  of 
Figo  (2d)  is  determined  by  transition  frequencies  attribut- 
able to  the  rotational  degrees  of  freedom  of  the  component 
gases.  Where  dispersion  of  this  type  has  been  detected,  the 
transition  frequency  vat  atmospheric  pressure)  has  been  high 
in  the  megacycle  region. 
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FIG  2 CORREGT'ON  FACTOR  FOR  THE  VELOCITY  IN  AIR  , AS  A FJNGTIQN 
OF  TEMPERATURE  AND  PRESSURE  . THE  CORRECTION  FACTOR  IS 
THE  FRACTIONAL  CORRECTION  THAT  MUST  BE  APPLIED  TO 
EXPERIMENTAL  MEASUREMENTS  TO  REDUCE  THEM  TO  ‘STANDARD 
CONDITIONS.  THE  FOUR  CHARTS  APPLY  TO  DIFFERENT  FREQUENCY 
REGIONS 
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